function [simulation_result] = FullWaveformInversionSimulation(anatomic_part_name, inversion_parameter)

%输入参数：
% AnatomicPartName:     tibia 或者 tibia_and_fibula
% InversionParameter:   velocity  或  density
%输出参数：
% simulation_result  ：与输入参数匹配的仿真结果

switch anatomic_part_name
    case "tibia"
        if inversion_parameter == "density"
            fid = fopen("../const_file/tibia_density.bmp", "r", "b");
            file_content = fread(fid);
            simulation_result = cast(file_content, "uint8");
            fclose(fid);
        elseif inversion_parameter == "velocity"
            fid = fopen("../const_file/tibia_velocity.bmp", "r");
            file_content = fread(fid);
            simulation_result = cast(file_content, "uint8");
            fclose(fid);            
        end
    case "tibia_and_fibula"
        if inversion_parameter == "density"
            fid = fopen("../const_file/tibia_and_fibula_density.bmp", "r");
            file_content = fread(fid);
            simulation_result = cast(file_content, "uint8");
            fclose(fid);               
        elseif inversion_parameter == "velocity"
            fid = fopen("../const_file/tibia_and_fibula_velocity.bmp", "r");
            file_content = fread(fid);
            simulation_result = cast(file_content, "uint8");
            fclose(fid);   
        end
     otherwise
        warning('Unexpected input parameters.')
end

end